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Abstract 

The maximal conserved phase gradient is introduced as an order parameter 

to characterize the transition from phase- to defect-turbulence in the complex 

Ginzburg-Landau equation. It has a finite value in the phase-turbulent regime 

and decreases to zero when the transition to defect-turbulence is approached. 

Solutions with a non-zero phase gradient are studied via a Lyapunov analysis. 

The degree of "chaoticity" decreases for increasing values of the phase gradient 

and finally leads to stable travelling wave solutions. A modified Kuramoto- 

Sivashinsky equation for the phase-dynamics is able to reproduce the main 

features of the stable waves and to explain their origin. 
PACS numbers: 47.27.Cn, 05.45. +b, 47.27.Eq 
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Spatially extended chaotic systems have been recently the subject of several theoretical 
and experimental investigations [1]. In spite of remarkable progresses, a still open problem 
concerns the possibility to describe nonequilibrium phase transitions with concepts borrowed 
from statistical mechanics [1-4]. The complex Ginzburg-Landau equation (CGLE) is one 
of the most appropriate model to study such transitions, because is universal [1,5] and ex- 
perimentally relevant [6]: the dynamics of extended systems undergoing a Hopf bifurcation 
from a stationary to an oscillatory state is described by the CGLE [7]. Moreover, several 
physical, chemical and biological phenomena can be well reproduced through the CGLE [1]. 
Even in one dimension the CGLE displays a variety of dynamical regimes and phase tran- 
sitions [2-4,7-9]. Among them much attention has been devoted to the transition between 
two different type of chaotic phases [2-4,7,8]: namely, the phase-turbulent (PT) and the 
defect-turbulent (DT) regimes. However, this transition has been mainly studied from the 
DT side [2-4]. 

This Letter will be focused on the PT regime and in particular on the introduction of an 
order parameter for the above mentioned transition. To be more specific, let me write the 
one-dimensional CGLE as 

A t = (1 + i Cl )A xx + A- (1 -ic 3 )\A\ 2 A (1) 

where the parameters c\ and c 3 are real positive numbers, and A(x,t) = p(x,t) exp[iip(x,t)] 
is a complex field of amplitude p and phase ip. Of particular interest in the (ci, c 3 )-plane 
is the so-called Benjamin- Feir line (BFL) defined by c 3 = 1/ci, which identifies the linear 
stability limit for the plane-wave solutions of (1) [5]. The PT regime is encountered just 
above the BFL (i.e. for c 3 > 1/ci) [7]. In this state the chaotic behaviour of the field is 
essentially ruled by the dynamics of the phase. Moreover, the amplitude is always bounded 
away from zero, accordingly; the average phase gradient v = j- J L dxd x tp(x,t) is conserved 
for periodic boundary conditions [7]. Another important line (L\) separates the PT state 
from a more chaotic phase: the DT regime. In this state, amplitude-dynamics becomes 
predominant over phase-dynamics [2,7,8]. In particular, large amplitude oscillations are 
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observed, occasionally driving p(x } t) to zero, in which case the phase is no longer well defined 
and v is not conserved. The vanishing of the amplitude determines a so-called space-time 
defect. The density of defects 8d is a good order-parameter to characterize the transition 
from the DT to the PT phase; in fact its value is > in the DT regime, while it vanishes 
approaching the PT phase [2,4]. Conversely, none of the other parameters introduced so 
far (e.g. phase and amplitude correlations lengths, or the Kaplan- Yorke dimension density) 
reveals a clear signature of the transition [3,4]. 

In the PT regime, a relevant parameter to characterize a state of the CGLE is the value 
of z/, which is a conserved quantity in the absence of defects. For each value of v < 1 there 
exists a plane-wave solution of the CGLE 



where fi = c 3 — (ci + c 3 )u 2 . Below the BFL, these solutions are stable against long- 
wavelength instabilities for v 2 < (1 — Cic 3 )/(2(1 + c^) + 1 — Cic 3 ) [5], while they turn out to 
be all unstable in the PT regime. However, it is reasonable to expect that different v values 
will characterize distinct classes of chaotic or quasi-periodic solutions also in the PT phase. 
Nevertheless, the majority of analyses reported in literature has been devoted to solutions 
with v = 0. 

In this Letter, I determine the maximal value vm of the conserved phase gradient in 
statistically stationary states, i.e. in the limit L — > oo and t — > oo. I shall argue that this 
can be used as an order parameter, with vm = in the DT regime, with vm > in the PT 
state, and with a smooth change at the transition line L\. Moreover, a Lyapunov analysis 
reveals that for increasing iz-values the solutions are less and less chaotic. This behaviour can 
be explained assuming that the phase-dynamics is ruled by a modified Kuramoto-Sivashinsky 
equation (MKSE) [7,10-12]. For reasonably high iz-values travelling pulse-train solutions are 
found to be stable. The main parameters characterizing such solutions can be derived from 
the above mentioned MKSE. 

The analysis here reported is limited to a parameter region where it is known that 
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the CGLE shows a continuous transition from the PT to the DT regime [2,4], namely 
C\ = 3.5 and c 3 > \jc\. The integration scheme here adopted is a time splitting code 
where the integration of the spatial derivatives is not performed, as usual, in Fourier space, 
but instead is computed in x-space through a convolution integral [13,14]. The adopted 
integration parameters are: spatial resolution dx = 0.5; integration time step dt = 0.05; 
number of grid points N = L/dx ranging from 2048 up to 8192, and periodic boundary 
conditions. 

In order to study solutions with a non zero value of z/, an initial state with v ^ has 
been prepared and its evolution followed. Due to the absence of defects in the PT regime, 
v should be conserved, at least for small initial values. I have numerically verified that v 
is indeed conserved for all values below an upper limit vm- For v > vm defects eventually 
arise leading, after a readjusting time, to a i/-value smaller than vm- This behaviour can 
be understood by observing that the minimal /9-value decreases continuously for increasing 
i/, until, for v > um } the amplitude can eventually vanish with the consequent emergence 
of defects. In order to evaluate the maximal phase gradient accurately, the conservation of 
each vm value reported in Fig. 1 has been checked considering from 50 to 70 different initial 
conditions. In particular, each configuration has been followed for a time t > 10, 000, after a 
reasonably long transient had been discarded. Moreover, for each examined c 3 -value, some 
trajectories have been followed for longer times, typically of the order of t = 150, 000. No 
dependence of vm on L has been observed, considering chain lengths from L = 1024 up to 
L = 4096. Fig. 1 shows a roughly linear decrease of vm for increasing values of c 3 . A linear 
fit of the data gives vm — — 0.74c 3 + 0.57. Assuming that the linear behaviour extends up 
to the L\ line, a critical value = 0.76 ± 0.03 is obtained. This result is in reasonable 
agreement with the values obtained from the defect density reported in literature [2,4], as 
well as with our own measurements of 8d, which give Cg = 0.755 ± 0.002 [13,15]. 

As a further characterization of these solutions, I have evaluated the maximal Lyapunov 
exponent Am for several values of c 3 and v. I noticed multistability for v > 0, i.e. several 
coexisting attractors for each v value, different attractors being characterized by different 
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Lyapunov exponents. However, a characteristic common to all these initial conditions con- 
sidered is that Am tends to decrease with increasing u, except for small fluctuations (see 
Fig. 2). In particular, for c 3 < 0.5 and for sufficiently high values of u, non-chaotic states are 
found. Two different kinds of non-chaotic states have been observed: the former (type a) is 
spatially periodic with spatial wavevector q = u } formed of identical "pulse-like" structure 
of length Lp = 2ir ju\ the latter (type (3) shows essentially periodic regions separated by 
domains in which p is constant and the phase decreases linearly (see Fig. 3). The selection 
of these patterns depends on the initial conditions. Numerically I found that all non-chaotic 
solutions are of the form: 



where = /?(£) exp[itp (^)] is in general complex. Amplitude and phase can be written 
as 



with £ = x — vt. 

For spatially periodic patterns, the elementary pulses of length Lp are stable solution 
of the CGLE for a short chain of length L = Lp and with a phase gradient up = 2ir/Lp } 
i.e. the minimal non-vanishing phase-gradient for such a periodic system. These travelling 
pulses originate through a bifurcation from the plane- wave solutions (2). Moreover, they 
are stabilized in short systems, because of the long-wavelength instability cutoff. A more 
extensive and detailed study of these short-chain solutions will be reported elsewhere. Here, 
I want just to point out that solutions with a conserved up are no longer observable below 
a minimal length L m i n . The corresponding phase gradient ujj = 2ir / L m i n is an upper bound 
for um-, as shown in Fig. 1. For increasing c 3 , the value um is better and better approximated 
by ujj. A naive explanation of this fact can be given by assuming that in the proximity of 
um, only stable solutions of type a are observed. For these solutions the minimal value of 
the spatial period is obviously L m i n and the maximal possible phase gradient ujj. However, 




(3) 




(4) 
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in the limit v — > um } for c 3 < 0.5 solutions of type a coexist with those of type /3, while 
for c 3 > 0.5 slightly chaotic solutions, formed by an array of elementary pulses of different 
lenghts L p > L m i n} are observed. All these facts imply that ujj > vm- 

As a final point, I would like to explain why Am decreases with increasing v. This 
can be done by recalling that just above BFL, the dynamics is essentially ruled by the 
phase behaviour, since the amplitude can be considered as a "slaved" variable of the phase 
[7,8,12]. A phase modulation <f>(x,t) = ip(x,t) — vx = ipo(£) -\-wt on the plane wave solution 
(2) satisfies a MKSE [11], 

j>-n + n&t + nPd 2 j + n { ?{d x <f>f + q 3 <9^ + n^d^ + n^d^j = o (5) 

where Oi = 2v(c 1 + c 3 ), 0^ = c lC3 -l, n {2) = (ci + c 3 ), 3 = 2i/ci(l + c§), 0^ = cj(l + c 2 3 )/2 
and 0^ = 2ci(l + c 2 ). In order to check the validity of Eq. (5) for the dynamics of the 
phase for the CGLE in the PT regime, I derive from Eq. (5) an expression for the velocity 
v and the frequency u of the travelling solutions (3). Following [11], I obtain 

u = fi - 0< 2) < (d s Vo) 2 > +0i 2) < (df V>o) 2 > (6) 

v = fii + (0^ 2) < (d^of > -n 3 < (d 2 ^o) 2 > +Ol 2) < («9 s Vo) 2 <9fV>o >)/ < (d^o) 2 > (7) 

where < • > is the average along the chain and over several consecutive realizations. The two 
expressions can easily be obtained by noticing that for a solution of type (4) the temporal 
derivative of the phase can be written as (j> = u — vd^tpo(^) } while d x (j)(x) = d^tpo(^). Substi- 
tuting the temporal derivative into Eq. (5) and averaging both sides of the equation leads 
to Eq. (6). To obtain Eq. (7) both sides of Eq. (5) should be multiplied by d^tpo(^) before 
averaging. Several quantities appearing in Eq. (5) have zero average due to the periodic 
boundary conditions. Therefore the final expressions for u and v are drastically simplified. 
By inserting the phase values obtained from simulations of the CGLE into Eqs. (6) and 
(7), a very good agreement with the measured quantities is indeed achieved (see Table 1). 
Nevertheless, I expect that the phase description (5) will become less and less accurate ap- 
proaching the L\ line and will finally break down when the amplitude dynamics becomes 
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predominant. However, I have verified that in the whole examined parameter interval the 
amplitude dynamics is essentially ruled by the phase, except when defects arise. 

Assuming that Eq. (5) describes sufficiently well the dynamics of the CGLE-phase, let 
me now explain the observed non-chaotic behaviours. In Ref. [II] it has been shown that 
Eq. (5), rewritten for the variable d x (f>, belongs to a class of equations that shows stable 
solutions formed by pulse-like periodic structures travelling along the chain with a finite 
velocity [16,17]. The birth of these non-chaotic solutions is due to the presence of the 
dispersion term oc d^(f> [16]. In particular, in Ref. [17] it has been shown that this equation 
is multistable: chaotic and non chaotic solutions can coexist. However, stationary periodic 
attractors prevail over the strange ones for increasing values of the dispersion constant 
fi 3 and, above a threshold value, only non chaotic solutions are observed. Since fi 3 is 
proportional to z/, this explains the decreasing behaviour of Am for increasing z/, as well as 
the origin of the observed travelling wave solutions for high v- values. Finally, it should be 
remarked that for c 3 > 0.5 it is no longer possible to observe stable solutions, because they 
would occur for value of v > vm- 

In conclusion, I emphasize that vm turns to be a good order parameter to describe the 
transition from the PT to the DT phase. The role of this new parameter is complementary 
to that of the defect density, since it is non-zero in the PT regime, while the defect density 
is non-zero in the DT state. Moreover, the dynamics in the PT phase depends strongly on 
the z/-value: for v — > a predominance of chaotic solutions is observed, while for v — > vm 
non-chaotic behaviours prevail and a new family of stationary solutions emerges. 
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FIGURES 



Fig.l: Maximal phase gradient vm as a function of the parameter c 3 (circle). The values 
for ujj are also shown (squares). The solid line represent a linear fit for the vm data. 

Fig. 2: Maximal Lyapunov exponents versus the phase gradient v for two different sets of 
initial conditions for c 3 = 0.5. The two initial conditions correspond to states with phase 
gradient v and with noise added only on the amplitude (circle) or added on both amplitude 
and phase (asterisks). The data have been obtained for system size L = 1024 and for 
integration times t = 130,000 - 250,000. 

Fig. 3: Phase and amplitude for the two kinds of observed non chaotic solutions: a) solution 
of type a with c 3 = 0.5 and v = 0.184; b) solution of type (3 with c 3 = 0.35 and v = 0.307. 
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TABLES 

TABLE I. The parameters characterizing the stable solutions (3), namely the propagation 
velocity v and the frequency u, are reported for various v values. The values are estimated using 
the expressions (6) and (7). In parentheses are the values measured directly from the travelling 
waves solutions. The data in the table refer to C3 = 0.5, except for (*) that is for C3 = 0.4. The 
label (+) indicates solutions corresponding to different initial conditions, but with the same v. 



V 


UJ 




V 


V 


UJ 




V 


0.123 


0.428 (0.427) 


1.26 


(1.24) 


0.123W 


0.428 (0.424) 


1.06 


(1.09) 


0.184 


0.355 (0.341) 


1.59 


(1.65) 


0.184(+) 


0.344 (0.338) 


1.58 


(1.63) 


0.245 


0.229 (0.223) 


1.93 


(1.98) 


0.282M 


0.074 (0.070) 


2.37 


(2.34) 



11 



